A quantum solution to Gibbs Paradox with few particles 
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We present a fully quantum solution to the Gibbs paradox (GP) with an illustration based on a gedanken 
experiment with two particles trapped in an infinite potential well. The well is divided into two cells by a solid 
wall, which could be removed for mixing the particles. For the initial thermal state with correct two-particle 
wavefunction according to their quantum statistics, the exact calculations shows the entropy changes are the 
same for boson, fermion and non-identical particles. With the observation that the initial unmixed state of 
identical particles in the conventional presentations actually is not of a thermal equilibrium, our analysis reveals 
the quantum origin of the paradox, and confirm the E. J. Jaynes' observation that entropy increase in Gibbs 
mixing is only due to the including more observables measuring the entropy. To further show up the subtle role 
of the quantum mechanism in the GP, we study the different finite size effect on the entropy change and shows 
the works performed in the mixing process are different for various types of particle. 



PACS numbers: 05.70.Ln, 05.3O.-d, 03.65.Ta, 03.67.-a 

Introduction- More than one century ago, J.Willard Gibbs 
[1] pointed out that, if the entropy were not extensive due to 
the neglection of indistinguishability of particles, there would 
be an entropy increase after two ideal gases from two con- 
tainers at thermal equilibrium are mixed. On the contrary, the 
inverse process would bring the system back to its original 
state and cause the entropy decrease for the closed system. 
This apparently violates the Second Law of Thermodynam- 
ics(SLT). However, for distinguisable particles, this inverse 
process would not bring the system to the original inital state, 
and not violate the SLT. 

In many textbooks, e.g., the Ref.[2], the GP is claimed to be 
solved by considering the indistinguishability of particle in the 
view of quantum theory, and adding correction to the expres- 
sion of entropy. This mixing process will bring no entropy 
change for identical particles, while additional increase for 
non-identical ones. However, this pedagogic explanation in- 
cludes neither the internal properties of particle nor a quantum 
peculiarity of the problem (e.g., non-commutativity). On this 
aspect, we still do not understand how the entropy changes 
disappear only when two gases are the same, even though they 
can be arbitrary similar. In the classical thermodynamic, it 
is always able to assumed the substances empirically distin- 
guished, however without defining how small their difference 
is. Thus, it is expected the perfect resolution of the GP should 
be referred to the quantum mechanics with intrinsic consid- 
eration of indistinguishability, even with internal variables of 
particles. This point was considered by some authors [3, 4] 
by introducing the internal state phenomenologically, and the 
discontinues of entropy change is avoided by utilizing another 
thermodynamic quantity- work. 

Another resolution was based on the recognition that the 
knowledge of difference in particle types serves as informa- 
tion, which enables extracting more work to compensate the 
entropy increase [5]. However how the internal state involve 
in the mixing process is still unknown. On this sense, a fully 
consistent quantum description of the GP, to our best knowl- 
edge, is still lack, especially on the indistinguishability of the 



particles based on the standard quantum theory, such as the 
second quantization approach with symmetrizing wave func- 
tions. Fortunately, the quantum thermodynamics (QT) [6-8] 
has enlighten solutions to some paradoxes in thermodynam- 
ics [9, 11], e,g, the Maxwell's demon (MD) paradox [10-12] 
with quantum Szilard engine. The QT may only concern the 
few-particle system, which sometimes is efficient to reveal the 
underlying principles and intrinsic mechanism. 

In this letter, we give a fully quantum solution to the GP 
involving the symmetrization of wave function of identical 
particles. We consider the mixing process of two particles 
confined in a infinite square potential well, which is divided 
into two sub-cells by a solid wall. Using the exact expres- 
sion of the density matrix, we demonstrate explicitly that the 
initial un-mixing systems in the conventional presentation of 
GP are not in thermal equilibrium with respect to the whole 
well. For the identical particles, we follow the standard ap- 
proach of second quantization and correctly write down the 
thermal equilibrium state of an un-mixing whole system by 
including internal variables and symmetrizing wave-function. 
This doing surprisingly restores the original version of GP 
for both non- and identical particles (Boltzman particles, and 
bosons/fermions) and solves the GP by showing the same en- 
tropy changes. The detailed comparisons with the case ignor- 
ing internal state is made to explicitly explain the origin of the 
paradox. For the cases with more particles, we also demon- 
strate almost the same result except for a little difference be- 
tween boson and fermion, due to the finite size effect. 

Gibbs paradox with infinite potential well- In thermal equi- 
librium, the the most probable state of many particles in a 
containers with two sub-cells of equal volume usually is the 
situation with each cell containing the equal number of parti- 
cles. Thus, to show up the essence of the GP, we firstly con- 
sider the two particle case, similar to the single molecule case 
[13] for the Szilard heat engine in resolving the MD paradox. 

Let us model the whole container as an infinite high square 
potential well V(x, I) (= for x E (0, 1) and oo other- 
where) to trap two particles (see Fig. 1(a) ). The well is 
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FIG. 1. (Color Online) Gibbs Paradox modeling with two particles 
trapped in a one-dimensional square potential well of infinite depth 
with width /. (a)The ordinary potential. The box is separated by 
a central ^-potential; (b) Different potentials trap different colored 
particles where the colors characterizing their internal states. The 
mixing is realized by isothermally moving the barrier of each poten- 
tial. 



initially separated by the infinite high ^-potential, namely 
vS (x — 1/2) with v — > oo. With the single particle eigen- 
wavefuctions \<Pn) ( on left (right) cell correspond- 

ing eigen-energies E% (E^) = 2h 2 n 2 n 2 /Ml 2 , the eigen- 
wavefuctions of two distinguishable particle system read as 
l^nro) = IVn)i ® \fm)2 corresponding to eigen-energies 
E^ v m = E% + E v m , for u,v = R,L. With the explicit ex- 
pressions of eigen-wavefuctions given in the Supplementary 
Materials (SM), (p^(x) ~ (x \ fn) and f^ix) are supported 
over the whole domain (—00,00) including (0,/). This ob- 
servation seems trivial, however can serves as the key point to 
resolve the GP as well as the MD paradox [11, 12] . 

To show this, we write down the thermal equilibrium state 
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of two Boltzman particles with inverse temperature (3 and par- 
tition function Zy. It includes not only the terms with 
and 1$^) for two particles in the different sub-cells shown 
in Fig. 1(a), but also the terms with \<&nm) and \$nm) f° r 
two particles in the same sub-cells. If the later terms were 
ignored as in the conventional description of GP, the initial 
un-mixing state used here would not be in thermal equilib- 
rium, and thus could not correctly address the very problem 
in GP correctly. For non-identical particles, we can actually 
assert that one specified atom in the right sub-cell and another 
in the left. As emphasized by Jaynes, the specified assertion 
to distinguishable particles actually reduces the system's en- 
tropy. Such difference between non- and identical particles is 
the very root of the GP. Based on this observation, we solve 
the GP by introducing internal state to the identical particles 
to restore the original GP in the next section. 

Quantum description - In order to resolve the GP, we would 
make the situation with equal number of particles in the both 
sub-cells be the thermal state by referring certain prior infor- 
mation. It is possible for non-identical particles since they 
could be subject to different potentials. With the similar 
method, we can prepare the initial thermal state for identi- 
cal particles with internal variable. Using the "colors" to label 
the internal degenerate states, the blue \b) and the red ones 



|?'). Let the blue and red particles be trapped respectively in 
the left and right sub-cells with the Hamiltonian 

H = p 2 /2M+\b) (b\ V(x, l/2) + \r) (r\ V(x-l/2, 1/2), (2) 

where the potential is shown in Fig. 1(b). And we assume that 
the internal state of the two identical particle are different. Ob- 
viously, this assumption assigns one bit information, which 
brings the identical particle with the same initial un-mixing 
state as non-identical one. Therefore, the Gibbs mixing could 
causes entropy increase for both non- and identical particles. 

The single particle eigen-wavefunctions of the present sys- 
tem are \bip^) = |&)<8)|<^) and \rip^) = Ir)®]^^) with cor- 
responding eigen-energy E^ and E^. Associating with their 
quantum statistics properties, the two-particle wavefunctions 
before the Gibbs mixing are symmetrized or anti-symmetrized 
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with the same eigen-energy En m Bi * F ^ — E% + E.^ 
Here, the sign plus/minus corresponds the boson/fermion 
case. These eigen-functions determine the density matrix 
of this system Pu-b(f) (f° r the explicit form see the SM), 
similar to Eq.(l) with the partition function Z U _ B ( F < ) = 

E„, m exp(-/^-^). 

The mixing process is carried out by moving the right side 
of left potential to the position x = I and the left side of 
the right potential to the position x = isothermally, illus- 
trated in Fig. 1(b). After the mixing process, two particles 
are constrained in the larger domain (0 < x < I). The single- 
particle wavefunctions of this system then become = 
\s, 4> n (x,l)} (degenerate for s — b,r,) with the same eigen- 
energy E„ = n 2 TT 2 h 2 /2M I 2 . The two-particle wavefunc- 
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tions after mixing 
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corresponds to eigne-energy En m B ^ = E n + E rn , which 
gives the density matrix Pm-b(f) with partition function 
Z M -b{F)- The entropy change in the mixing process is sim- 
ply calculated by using von Neuman entropy as ASb(f) — 

Tr [Pu-B(F) InpM-B(F)] - Tr [p M -B(F) In Pm-b(f)]- 

We notice that the calculations about entropy change con- 
cern only the eigen-energies, rather than the concrete form of 
the wavefunction. However, the symmetry type of the wave- 
fucntion determines the counts of states to the partition func- 
tion. For two non-identical particles with un-symmetrized 
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eigen-wavefunctions \^n m N ) = \b(Pn)i®\ rl Pm 
sponding eigen- energy simply is E^ N = E% - 
using the wavefunctions after mixing hl , *^" iv ) 
\r4>m)2 w tth the eigen-energies E^f^ — E n 
tain the same expression of the entropy increase AS^. There- 
fore, it is concluded that the entropy changes in the process of 
mixing two bosonic, Fermionic and non-identical particle are 
the same, namely 
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FIG. 2. (Color online) Entropy change in the process of mixing two 
atoms as a function of (a) trap size I with different inverse tempera- 
ture f3 = 0.5 and (b) inverse temperature j3 with different I = 10 for 
different species: Fermi (Blue-solid), Bose(Red-dashed) and Non- 
identical(Green-dashdotted); Work done as function of (c) trap size 
and inverse temperature(d). The classical case of entropy change 
ASciassicai is also marked as gray-dashed line in subfigure (a) and 
(b). 



On this sense, GP is resolved. 

It is also meaningful to compare the present results with 
that without considering the internal state for Bosonic and 
Fermionic particles. In our consideration, the initial un- 
mixing state is a thermal state with all the three situations 
in Eq.(l). The mixing process is completed by removing the 
central potential isothermally. The entropy change AS' X and 
work W' x — ksT [lnZ^_ x — lnZ' AI _ x ~\ are calculated by 
using the partition functions of un-mixing and mixed parti- 
cles Z' u _ x and Z' M _ X with X — B,F,N. Here, we deal 
with the usual setup in the conversional presentation of GP 
for non-identical atoms, which is the same as that with inter- 
nal state. We refer the detailed analysis for different species 
to Supplement Materials. One prominent feature is that both 
entropy changes and work done for different species are dif- 
ferent, while these are the same in the process above with the 
correct considerations about internal state. We show these dif- 
ferences in Fig. 2. In Fig. 2(a-b), the entropy changes of 
identical particle approaches zero, while that of non-identical 
particles is 2 In 2. The discrepancy just recovers the discon- 
tinuousness of entropy changes in conventional presentation 
of GP. The similar difference in work extraction is illustrated 
in Fig.2(c-d). In the present model, no internal freedom is 
probed and no prior information is added for identical par- 
ticles. In comparison with the previous presentation of GP 
with internal freedom, we can conclude that the discontinu- 
ousness is caused by the different start-point between identical 
and non-identical particles. 

To understand the above result, we recall the conventional 
presentation of GP that the initial state of the identical parti- 
cles is formally different from the non-identical ones, which 
results in no entropy change for mixing identical particles. 
To our understanding, the same entropy changes are restored, 



solely by restoring the GP problem with the same start points 
for the different types of particles. With this observation, we 
conclude that the GP is rooted in the initial difference of un- 
mixing state rather than the character of indistinguishability 
of particles. It is also noticeable that we only utilize the quan- 
tum definition of entropy, the von Neuman entropy, other than 
any phenomenological presentations. 

Mixing entropy and work of two particle systems - As for 
the quantum effect in our approach for GP, the finite size of 
the well induces many interesting phenomenon. To explicitly 
demonstrate the effect, we take the bosonic case for exam- 
ple. In terms of the Theta function 63 (0, q) = 1 + 2 ^ ra 1 n < 
the partition functions before and after the Gibbs mixing are 
expressed[14] as Z V -b = [O3 (0,q) ~ l] 2 /4 and Z M -b = 
[6*3 (0, q 1 / 4 ) — l] /4 respectively with the parameters q — 
exp (— 2f3n 2 h 2 /ml 2 ) . The straightforward calculation ex- 
plicitly gives the entropy increase during the mixing process 
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of the Theta-fucntion 

f? 3 (0,exp (7r 2 /lng)) ,we 



Using the duality 

(-lng/7r) 1/2 (93(0,g) 

prove that #3(0,(7) — > l/v~ m l ?/ 7r m tne high tem- 
perature limit T — s- 00 or in the classical limit L 00. 
With this observation, the classical result is recovered as 
linii-s-oo ASs = 2 In 2 = AS c i asS i ca i. 

Very quantum nature is the dependence of entropy change 
on size of the trap and also on the temperature, illustrated 
in Fig. 3(a). The entropy change A5 c iassicai in the classi- 
cal case is marked as gray-dashed line in Fig. 3(a-b). The 
entropy changes tends to the classical one AiS G i ass i ca i as the 
trap size approaches infinite, presented in Fig. 3(a). This con- 
firms our theoretical analysis given above. We also show its 
dependence on the inverse temperature /3 in Fig. 3 (b). The 
entropy also tends to the classical case as the inverse temper- 
ature approaches zero(/3 — > 0), since the thermal fluctuation 
compares well the discreteness of the energy levels in the well. 
Finally, starting with the correct uses of the eigen- vectors and 
-energies, we consider work done in the mixing process as 
the references[4, 5]. The mixing process is performed isother- 
mally, and the work just compensates the free energy change, 
namely, Wb = ksT [In Zjj_b — In %m-b\ which is the 
same for different species. This seemingly-trivial observation 
also solves GP. We show in Fig. 3(c-d) the work done as func- 
tion of trap size I and inverse temperature /3 respectively. With 
larger trap size, the work approaches zero. However, it is not 
zero as temperature increase. Similar behavior has also been 
observed in the insertion process of SHE[11]. Theoretically, 
the work diverges as \JT as T — > 00. The detailed derivation 
is presented in Supplement Materials. 

Mixing entropy of 4 particle systems- To give prominence 
to the properties of quantum statistics, we consider the mix- 
ing process of four particles. Initially, there are two particles 
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FIG. 3. (Color online) Entropy change in the process of mixing two 
atoms as a function of (a) trap size I with different inverse temper- 
ature ft = 1 (blue-solid), 2 (red-dashed) and 3 (green-dashdotted) 
and (b) inverse temperature /? with different I = 10 (blue-solid), 20 
(red-dashed) and 30 (green-dashdotted). Work done as function of 
(c) trap size and inverse temperature. The classical case of entropy 
change is marked as gray-dashed line in subfigure (a) and (b). 



in each sub-cell. The internal states of the particles in the 
right/left side are all blue/red. Its partition function in thermal 
equilibrium is initially [14] 



Zu-B(F) (N = 4) 



Z, 



U-B{F) 
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where Z" ( } {q) = ±[(Zi(q 4 )) ± Z^q 8 )} denotes the par- 
tition function of a two-particle system confined in a square 
potential well with width 1/2. For Z x {q) = [0 3 (0, q) - 1] /2, 
there exist two limit results: in the quantum limit with 
low temperature j3/l 2 — > oo, we have Z\{q) ~ 2q; 
and in the classical limit f3/l 2 — > 0, we have Z\{q) ~ 
1/2 \J — \nq/n— 1/2. Together with the partition function[ 14] 

Zm-b(f) (N = 4) = (zf (F) (<7 1/4 )) 2 after mixing, we cal- 
culate the entropy change is 
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In the classical limit, we have lirm_ i . 00 ASb(f) (N = 4) = 
4 In 2. Indeed, when deviating from the classical limit or the 
lower temperature, the entropy change obviously displays the 
obvious differences between boson and fermion. 

If we mix identical particles without considering 
the internal state, one must start from an equilibrium 
state with partition function[14] Ztj_b(f) (N = 4) = 
En=o Zn (F) (q)Z^(q), and that Z M _ B{F) (N = 4) = 

ZP \q 1/4 ) after mixing. In the classical limit, we obtain 
the entropy change as lirri^oo AS B , F s (N = 4) = 0. 
The results obtained here can be directly generated for the 
arbitrary atom number N, 



Summary-We have resolved the Gibbs paradox in a fully 
quantum framework with the correct presentations of the ini- 
tial thermal states of the working substances consisting of the 
particles with different quantum statistical properties. Here, 
we utilize a standard quantum description, the wave function 
symmetrization from the second quantization of the particles 
with internal variable. The key point is our finding that the 
problem in the paradox are rooted in the different uses ( some- 
how misuses) of the initial thermal state for identical and non- 
identical particles. We use examples with two and four parti- 
cles to illustrate our comprehensive understanding for GP. We 
show that the entropy change of identical and non-identical 
particles are the same for two particle system, but the defer- 
ence could only be found as the finite size effect in the low 
temperature, or in the cases with more than two particles. 
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Supplementary Material 

In this supplemtary material, we provide the detailed calculations of the Letter. 

I. Wavefunctions 

The un-mixing single particle wave-function is listed as 

<^> = l^ Sin(2nW0 Q<X<1 '\ (7) 
\ otherwise 

/ , r\ f \/|sin[2n7r(a;-Z/2) II] l/2<x<l 

\AVn) = \ ^ • (8) 

[ otherwise 

with the same eigen-energy E% = = 2n 2 ir 2 h 2 /Ml 2 . After the removing the central potential, the atoms can stay in the 
whole area. The wavefunction is 

<^} = {V^ Sin(nW0 °< x < 1 , (9) 
[ otherwise 

with energy E n = n 2 ir 2 h 2 / (2M1 2 ). 

II. Mixing process of 2 atoms 

(i) With internal Freedom 

The Hamiltonian for the atoms with internal freedom is written as 

H =\b){b\®H b + \r){r\®H ri (10) 

where Hi = p 2 /2m + Vi(i = b, r) is the corresponding Hamiltonian of the atom with internal state i = b,r. For this model, the 
mixing is done by moving the right boundary of the left potential to the right (x = I) and the left boundary of the right potential 
to the left (x = 0). For the Bosonic case, the partion function of unmixing two atoms is 



Zu_ B = J2 cxp [-2 f3n 2 h 2 (n 2 + m 2 ) /Ml 2 



n.m—1 

\„=1 / \m=l / 

=Z 1 (q A ) 2 , (11) 

where q = exp (-Pir 2 H 2 /2Ml 2 ) and Z x (t) = £~ =1 r" is the one-particle partition function and it can be related to the Theta 
function by 

Zi(t) = [0 3 (O,t)-1]/2. 
The partition function of mixed atoms is obtained similarly as 

Z M -B = Z 1 (q) 2 . (12) 

The results for Fermionic and non-identical atoms is the same as that of Bosonic case. It is clear that q e (0, 1) and approaches 
1 as T — >• oo or / — >• oo. The Theta-fucntion has a duality feature 

V- (In q) Jw 9 3 (0, q) = 9 3 (0, exp (ir 2 / In q) ) . (13) 

With this equation, we explicitely get 

i \ -V 2 



(14) 



l \ - 1 / 2 



1 + 2 V] exp 



n=l 
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The factor exp (71" 2 / In q) is no more than 1 when T or lis sufficiently large, which leads to 

n=l v ' n=l v 7 

Thus, Y^=i ex P (inf) a PP roac hes as T — > 00 or Z — > 00 and the aympotic behavior of 3 (0, q) as T — > 00 or I — > 00 is 

3 (O,9)= f- 1 ^ 7 [l + o(l)]. (16) 



7T / 

Using the above relationship, the classical limit can be obtained as 



( I M l 2 1 



lim A 5b = lim 

Z^oo L^oo V' 013 I ( l 2 Ml 2 ,V 

= 2 In 2. (17) 

(ii) Without internal freedom 

In this section, we will discuss the mixing process of two identical atoms without any internal freedom. In this suituation, the 
two atoms has the probability to stay in the same side of the chamber, illustrated in Fig. l(e-f). 

Bosonic Case: The eigen-functions before mixing are \&nm B ) RR = [\ ( Pn) 1 |</>m) 2 + \'Pm) 1 Wn) 2 ] \®n^ B ) RL = 

[Wn) 1 \<Pm) 2 + Wn) 2 ] I ^ aIld \®n m B ) LL = [Wn) 1 |^m) 2 + \<Pm) l K) 2 ] fovn^m and when n = TO the 

eigen-functions |6^ B ) RR = \ ( Pn) 1 |¥>«) 2 an ^ = l^n)i Wn)^ The partition function is then derived as 

00 

Z 'u-b= E (2 + S nm )exp[-2l37r 2 h 2 (n 2 + m 2 )/(Ml 2 )] 

n.m—1 

^l&^q^-lf + l&^q*)-!]. (18) 

After the mixing, the two-particle wavefunction are jO^f^ 5 ) = [|0n)i l<?W 2 + 10m)! I</Vi) 2 ] / \/2 for n 7^ to and when n = m 
\®nn B ) = l^n)i I0n) 2 - The partition function is 

00 

^m-b= E 2 (l + Ucx P [> 2 fi 2 (n 2 +m 2 )/(2Mi 2 )] 

n,m=l 

= 1 -[9 3 (0,q)-l] 2 + 1 -[9 3 (0 1 q 2 )-l]. (19) 

Fermionic Case: The function before mixing should be anti-symmeterized as \®nm F ) RR = 
M)i K) 2 - K>i K) 2 ] l -n F ) iL = K) 2 - K) x K) 2 ] /V2 with n ^ to and \@ u n ~ F ) RL = 

M)M\-\^)M)^ /vs. 

Z' V _ F = ( 2 - <W) exp [-2/?^ 2 ft 2 (n 2 + to 2 ) / (Ml 2 )] 

nm 

= \[0s(0,<l 4 )-l] 2 -\[03(0,<l 8 )-l]- (20) 

After the mixing, the two-particle wavefunction are \<dnm F ) = [\<j> n )i \<fim) 2 ~~ \4>m)i |<An) 2 ] with n 7^ to. The partition 
function is 

Z M-F = E ^ ( X - ^ ° X P t-^ 2fi2 + / ( 2mL ')] 

nm 

= 1 -{8 3 (0,q)-l] 2 - 1 -[e 3 (Q,q 2 )-l]. (21) 



III. Mixing process of 4 atoms 
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n = 2 


i (Z 1 (q) 2 ±Z 1 (q 2 )) 


n — 3 


1 (Zi(g) 3 ±3Zi(g 2 )Zi( 9 ) + 2Z 1 (g 3 )) 


n = 4 


! / Z 1 (g) 4 + 3Z 1 ( g 2 ) 2 ±6Z 1 (g 4 ) \ 
24 ^ ieZifo 2 )^) 2 + 8Z 1 (q s )Z 1 (q) J 



TABLE I. The relationship between Z n (q) and Z\ (q) for n = 2,3,4. 



(i) With internal freedom 

At beginning, there are two blue atoms (whose internal states are \b}) in the left compartment and two red atoms (whose 
internal states are \r)) in the right compartment. Thus, the initial partition function is 



.. B(F) (JV = 4) = (z^ F \q))\ (22) 



where Z^ F ^ (q) denotes the connonical partition function of a two-particle system confined in a square protential well with 
width 1/2. The explicit expression of Z^ F \q) is 

Z^ F \q) = \{{Z x {q')f±Z l ^)). (23) 
The concrete form of the un-mixing gas is rewritten as 

Zu-B(F) (N = 4) = i ((ZiG? 4 )) 2 ± Z^q 8 )) 2 



I //fy, ( o,^)-iV ± g 3 (0,g 8 )-i \ (24) 



4 

After mixing, these four atoms are confined in a infinity high square well with width I, The partition function is simply 

2 



Z M _ B{F) [N = A) = (z^ F \q^))' 



4 

1 ( fe 3 (0,q)-l\\ 9 3 (0,q 2 )-l 



2 

, . . , , ± ^ 1 ! • (25) 

Therefore the entropy change is 

ac ci aa m Z m-b(f) (N = 4) 

AS* s(i?) = (1 - /S^) In— — — . (26) 

Zu-B(F) — 4) 

Using the property of Theta function, it is not difficult to show that under the classical limit (T — > oo or I — > oo ), the entropy 
change tends to be 4 In 2 for both Boson and Fermion systems, 
(ii) Without internal freedom 

In this section, we discuss the mixing process which is the same as the one in the last section except that the atoms considered 
here do not contain internal freedom. In this circumstance, the initial state is a mixing of 5 situations, i.e., the n-th situation is 
that there are n atoms in the left compartment while the others are in the right compartment. Thus, the initial partition function 
is 

4 

Z'u-BiF) {N = 4) = Z^ F \q)Z i _ n {q), (27) 



n=0 



where Z„^ F> (q) denotes the connonical partition function of a n -particle system confined in a square protential well with width 
1/2. The relationship between Zn^ F \q) and Zi{q) can be found in Tab.(I) . 
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After the mixing process, the system becoms to be a well containing four identity atoms and the partition function is 

7' („W\- 1 ( Z 1 (Q)+3Z 1 (q 2 ) 2 ±6Z 1 ( q i ) \ 

Zu - B ^ (N - A) - Uq ] -^[ ±6Z l{q >)Z l{q r + 8Z l( ,3)Z l(9 ) ) • (28) 

This partion function can also be related to the Theta function by the ways mentioned above. The entropy change during this 
process is 



z 'm mf\{ n = 4) 
AS B{F) (N = 4) = (1 — Wp) In f" B(F) 

U—B(F) ^ V — 4 J 



which can also be proved that lirn^oo AS B , F -. (N = 4) = by using the property of the Theta function. 



